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We study the consequences of deterministic chaos for diffusion-controlled reaction. As an example, 
we analyze a diffusive-reactive deterministic multibaker and a parameter-dependent variation of it. 
We construct the diffusive and the reactive modes of the models as eigenstates of the Frobenius- 
Perron operator. The associated eigenvalues provide the dispersion relations of diffusion and reaction 
and, hence, they determine the reaction rate. For the simplest model we show explicitly that 
the reaction rate behaves as phenomenologically expected for one-dimensional diffusion-controlled 
reaction. Under parametric variation, we find that both the diffusion coefficient and the reaction 
rate have fractal-like dependences on the system parameter. 

PACS numbers: 05.45. +b, 05.60,+w, 47.53. +n, 47.52. +j, 47.70.-n, 82.20.-w, 82.30.-b 



Matter is most often the stage of reactions which evolve on a reactive time scale being intermedi- 
ate between the long time scale of hydrodynamic transport phenomena and the short time scale of 
\ microscopic chaos. This chaos is generated by the collisions between the atoms and molecules of the 
■ fluid beyond a temporal horizon caused by the Lyapunov instability of motion. Under nonequilibrium 
0^ ' conditions, long-time trajectories organize themselves in phase space to form fractal structures and 
weird invariant or conditionally invariant measures, the consequences of which have just started to be 
explored. In this perspective, we study here simple models of diffusion-reaction processes in order to 
confront the phenomenology with the new approach based on deterministic chaos. 



I. INTRODUCTION 



>V „ , 

Irreversible phcnomenological equations such as the Navier-Stokes or diffusion-reaction equations describe transport 

and reaction processes in fluid flows or chemical reactions occurring on macroscopic spatio-temporal scales. For 
instance, spatial inhomogeneities of size C in the density are damped by diffusion over a time scale of the order of 
Tfjiff ~ C 2 /D, where D is the diffusion coefficient. In typical laboratory experiments, the size C of the inhomogeneities 
. is of the order of millimeters, centimeters, or larger so that the relaxation time Tdiff is macroscopic. 

Recent works have studied the relationship between these macroscopic processes and the much faster process of 
chaos in the microscopic motion of atoms or molecules in the fluid Jl|-0| ■ The defocusing character of the collisions 
between the particles of the fluid is at the origin of a very high dynamical instability in the microscopic motion. 
This instability, which is a major phenomenon at the center of the current interest, is characterized by a spectrum 
of positive Lyapunov exponents, as shown by many recent numerical and analytical studies |]§|-|icill- The maximum 
Lyapunov exponent of a dilute gas of particles of diameter d is typically of the order of 

v I 

Amax ~ -T In - , (1) 

t d 



where I is the mean free path between the collisions, and v is the mean velocity of the particles |11| , |12| . Consequently, 
the time scale over which the dynamical instability develops in the motion of a particle due to the collisions with 
the surrounding particles is of the order of the inverse of the maximum Lyapunov exponent which takes the value 
T c haos ~ ^max ~ 10~ 10 sec for a gas at room temperature and pressure. 

The microscopic time scale is in contrast with the macroscopic time scale of transport and reaction processes. This 
presence of two different time scales is an essential feature in the recent establishment of quantitative relationships be- 
tween macroscopic transport and microscopic chaos in both the thermostatted-system and the escape-rate approaches 
[p|^,P|Jl^dT5t| . Indeed, in both approaches transport coefficients are related to differences between two characteristic 
quantities of chaos: The thermostatted-system method works with the difference between the maximum and the 
absolute value of the minimum Lyapunov exponent ||,|],^]], whereas the escape-rate approach employs the difference 
between the positive Lyapunov exponents and the Kolmogorov-Sinai entropy pp|,p^Hj5[|- In this sense, the nonequi- 
librium property of transport is related to a slight disbalance between the dynamical instability, which is the cause, 
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and the induced temporal disorder, which is the effect. At equilibrium, the effect exactly compensates the cause. 
Away from equilibrium, temporal disorder is slightly reduced to the benefit of transport, which appears as an extra 
effect of the dynamical instability beside the temporal disorder. In this scheme, the fluid appears at the stage of a 
highly chaotic motion of its constituent particles, which animates different possible transport processes if the system 
is maintained out of equilibrium. 

The hypothesis of microscopic chaos [[12) or the chaotic hypothesis replace the old stochastic hypotheses in 
nonequilibrium statistical mechanics. Previously, the stochastic hypotheses assumed that transport processes arise 
from stochastic effects, such as the Langevin white noise which has an infinite Kolmogorov- Sinai entropy per unit 
time. However, these processes have thus a much larger temporal disorder than allowed by Newton's deterministic 
equations of motion [17]] . The new chaotic hypothesis has the enormous advantage of assuming a temporal disorder 
which is now compatible with the determinism of a microscopic Newtonian dynamics. In this regard, the chaotic 
hypothesis overwhelms the previous stochastic hypotheses, which nevertheless remains of great usefulness in their 
domain of validity. 

The purpose of the present paper is to describe several consequences of the hypothesis of microscopic chaos in 
diffusion-reaction systems. This class of physico-chemical processes has not yet been explored in the perspective of 
understanding their kinetics on the sole assumption of microscopic chaos, without involving stochastic Langevin or 
birth-and-death processes. The diffusion-reaction systems are particularly important in various fields of chemical 
physics, such as chemical kinetics, homogeneous and heterogeneous catalyses, pattern formation in nonequilibrium 
reactions and in morphogenesis, recombination processes in solid or liquid phases, as well as high-energy reaction 
processes in astrophysical systems Q . 

We shall focus here on the simplest diffusion-reaction process with a linear chemical reaction law, 

A ^ B , (2) 

which already provides a nontrivial dynamics. As a vehicle of our study we shall use multibaker models, which are 
deterministic versions of discrete Markov processes. The deterministic dynamics of a multibaker has therefore the 
same finite and positive Kolmogorov-Sinai entropy per unit time as the corresponding discrete Markov process. The 
models we propose are spatially extended generalizations of a baker-type model of isomerization previously studied 
by Elskens and Kapral ]l9| |. 

This paper is organized in the following way: In Section II, we construct multibaker models of diffusion-reaction by 
starting from a diffusive-reactive Lorentz gas. In Section III, we focus on the diffusive properties of the multibakers. 
By employing quasiperiodic boundary conditions we show that for the simplest model the diffusive properties are the 
same as those of the previously discussed dyadic multibaker of diffusion. We then demonstrate that under parametric 
variation of this model the diffusion coefficient exhibits a self-similar structure reminiscent of fractal curves. In Section 
IV, we describe the reactive properties. With quasiperiodic boundary conditions, we first study the simplest model 
and derive the dispersion relation of the chemiodynamic modes. We explicitly construct the phase-space distribution 
of these modes and define the reaction rate by comparison with phenomenology. We then show that the reaction 
rate behaves in a highly irregular manner if we consider the parameter-dependent model. Conclusions are drawn in 
Section V. 



II. DETERMINISTIC MODELS OF DIFFUSION-REACTION PROCESSES 



A. Definition of the models 



In order to motivate the introduction of the multibaker models, we first consider a reactive Lorentz gas in which a 
point particle undergoes elastic collisions on hard disks which are fixed in the plane. The disks may form a random or 
a regular configuration (Fig. |l|). A fraction of the disks are supposed to be catalysts where the point particle changes 
its state, or color, from A to B or vice versa at the instant of the collision. The mass of the particle is assumed to be 
the same in both states A and B. The phase space coordinates of each particle are given by its position, its velocity, 
and its color (x,y,v x ,v y ,c) with c G {A, B}. Since energy is conserved at the elastic collisions, the magnitude of 

the velocity is a constant of motion, v = \Jv x + Vy, so that the coordinates reduce to (x,y,tp,c) — (X, c), where 

ip = axcta,T\(v v /v x ) is the angle between the velocity and the x-axis. 

The motion induces a time evolution of the phase-space probability densities, or concentrations, for each color, 

f W = (/(x,f)) ■ ( 3 ) 
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The mean phase-space density, defined by the average of the concentrations, does not distinguish between the colors. 
Thus, we may expect that 

/(X) = ±[/(X,A)+/(X,i?)] (4) 



evolves in time exactly as in the non-reactive Lorentz gas, as studied elsewhere p0|-f22[. The dynamics of reaction 
should appear in the difference between the concentrations 

g(X) = f(X,A)-f(X,B) (5) 

which is expected to follow a macroscopic relaxation toward zero if there is equipartition of particles between both 
colors. 

It has been explained elsewhere that the flow dynamics of the Lorentz gas can be reduced to a Birkhoff mapping 
from collision to collision p2[ . Each collision can be represented by two variables: the angle 9 giving the position of 
impact on the perimeter of the disk as (x = cos 9, y = sin#) and the angle <p between the velocity after the collision 
and the normal at impact. The sine of the velocity angle 137 = sin <p together with the position angle 9 are the so-called 
Birkhoff coordinates, in which the mapping is area-preserving. All the collision events with the disk of label I are thus 
represented by the rectangle 

Tli = {(9, w, I) : < 9 < 2tt , -1 < w < +l} . (6) 
The dynamics of collisions can therefore be written as the Birkhoff mapping 

(0n+l, ln+1, Cn+l) = ni ^ ni^ni C-n) i (7) 

which is known to be area-preserving, time-reversal symmetric, and of hyperbolic character. 

A caricature of this mapping is provided by a multibaker model (Fig. ^) in which we suppose that the rectangular 
domains Eq. (|^) representing the disks are replaced by squares 

Si = {(x,y,l) : 0<x<l , < y < l} , (8) 

where I € Z is the label of the square. Each square of the multibaker model corresponds to a disk of the Lorentz gas. 
Now, the collision dynamics is simplified by replacing the complicated Birkhoff map Eq. (Q) by a baker-type map 
with horizontal stretching by a factor of two, followed by cutting the elongated square into two. The collisions from 
disk to disk are replaced by jumps of the particle from square to square according to the transition rule I — ► I — 1 if 
x < 1/2 and I — ► I + 1 if a; > 1/2 between next- neighboring squares. The squares are arranged such that they form a 
one-dimensional chain. One out of L squares is assumed to be a catalyst where the color changes from c = A (resp. 
B) to its complement c = B (resp. A). The map of the model is thus 

(2x,!,Z- l,tA , 0<a;< \ , 

2x - l,2±i,Z + l,c'J , | <x < 1 , 

where c' = c if I = 0, ±L, ±2L, ... and c' — c otherwise. This map is area-preserving, time-reversal symmetric, and 
chaotic with a positive Lyapunov exponent A+ = In 2 and a negative one A_ = — In 2, such as the multibaker map 

EE- 

In our model, the reaction is controlled by the diffusion if the reactive sites are diluted in the system. This important 
case of chemical reactions has been much studied in the literature since Smoluchowski seminal work p7[ . It is known 
that the macroscopic reaction rate is determined by the time taken by particles to diffuse toward the reactive site. 
A crossover occurs at dimension two which is the Hausdorff dimension of a Brownian path. Accordingly, the flux of 
reactants toward a catalyst is sensitive to the presence of the next-neighboring catalysts in less than two dimensions, 
but not in systems of dimensions higher than two. In particular, in a one-dimensional system the reaction rate should 
behave as 

D , , 

« ~ £2 ' ( 10 ) 

where D is the diffusion coefficient and L is the distance between the reactive sites or catalysts. 

A main goal of our work is to investigate the dynamical properties of our diffusion-reaction model in order to 
know whether this expected macroscopic behavior is confirmed from the microscopic dynamics or not. We shall also 
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consider a parametric variation of this model with a more complicated dynamics. This is due to an extra dependency 
on a shift parameter which is introduced when the half squares are glued back into the chain. When this continuous 
parameter varies it induces topological changes in the trajectory dynamics which are reminiscent of the topological 
changes induced by varying the disk radius in the Lorentz gas |20|-|22| . The parametric extension of the multibaker 
model shows that the diffusion coefficient as well as the reaction rate may vary in a highly irregular fashion as a 
function of a parameter. This is an important consequence of deterministic chaos which already appears on the level 
of one-dimensional maps, as will be discussed in the following section. 



B. The Frobeni us- Perron operator and quasiperiodic boundary conditions 



Since Boltzmann's work, it is well known that transport and reaction-rate processes should be conceived in a 
statistical sense because the individual trajectories are affected by the famous Poincare recurrences. We therefore 
consider the time evolution of statistical ensembles of trajectories as represented by the probability densities Eq. (||). 
They evolve in time according to the Frobenius-Perron equation 



ft+i(x,y,l,c) = ft * 1 (x,y,l,c) = (Pf t )(x,y,l,c) (t G Z) 



(11) 



We choose quasiperiodic boundary conditions by assuming that the solution of the Frobenius-Perron equation is 
quasiperiodic on the chain with a wavenumber k. Moreover, we suppose that the solution decays exponentially with 
a decaying factor \ = exps where \\\ < 1 or Re s < 0, 



ft(x,y,l,c) ~ x* exp(ifcZ) 



(12) 



The decay rate s is calculated by solving the eigenvalue problem of the Frobenius-Perron operator. We note that 
the Frobenius-Perron operator is in general non-unitary so that root vectors associated with possible Jordan-block 
structures may exist beyond the eigenvectors. We shall focus here on the eigenvectors because they control the slowest 
decay on the longest time scales p4| , p5| . 

For quasiperiodic solutions, the Frobenius-Perron operator reduces to the following Frobenius-Perron operator Qk 
which depends on the wavenumber k and acts on functions which are defined only in L successive squares of the chain, 



ft+i(x,y,0,c) = ( 

f t+ i(x,y,l,c) = ( 
ft+i(x,y,2,c) = ( 

ft+i(x,y,L- 2,c) 
ft+i(x,y,L- l,c) 



\-y) / t (§,2y,l,c) + e- ikL e(y-l) / t (£±I,2y - 1,£ 
\-y) / t (f,2y,2,c) + o(y - ±) f t (2^1, 2y - 1, 0, c) , 
\-v) /t(f,2y,3,c) + 6(y-i) / t (2±l, 2y - 1, 1, c) , 



l,c 



(13) 



9 {\-v) /t(f,2j/,i-l,c) + 6(y-\) /*(2±l,2y-l,Zi-3,c) , 
e lkL e(\-y) /t(f,2y,0,c) + 6{y-\) f t (^, 2y - 1, L - 2, c) . 



We notice that there is a reaction, i.e., a change of color, for particles passing the cell / = so that a concentration 
with I = and c appears in the second and in the last line. We also notice that there is a factor exp(— ikL) in the 
first line for the particle coming from the previous segment of length L in the infinite chain, where the concentration 
functions are multiplied by the factor exp(— ikL). On the other hand, there is a factor exp(ikL) in the last line for 
the particle coming from the next segment where the concentration functions are multiplied by exp(zfcL). Otherwise, 
this Frobenius-Perron operator is the same as in the infinite dyadic multibaker model studied in Refs. |23|-|25|l . 

As we discussed in the previous section, the presence of two chemical components c = A or B implies that the 
Frobenius-Perron operator Qk acts on 2L functions which can be linearly combined to separate the functional space in 
two subspaces on which two decoupled Frobenius-Perron operators would act. The first subspace is defined by Eq. (Q) 
where the Frobenius-Perron operator reduces to the diffusive Frobenius-Perron operator of the multibaker map. The 
second subspace is defined by Eq. (||) which gives a different evolution operator of reactive type. Diffusive properties 
are studied in the next Section III, while reactive properties will be discussed in Section IV. 



III. DIFFUSION DYNAMICS 
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A. Diffusive modes of the dyadic multibaker 



In this subsection, we consider the diffusive dynamics of the dyadic multibaker model Eq. (|9|) with quasiperiodic 
boundary conditions. The subspace of diffusion is defined by the mean density of Eq. (^). For hydrodynamic modes 
of wavenumber k we thus write 



f t (x,y,l) = ^[ft{x,y,l,A) + f t {x,y,l,B)] = exp 



and the new function obeys the simpler evolution equation 



i[ k + 2n- ]l 



Vt+i 



Vt(x,y) 



1 



(14) 



-i(k+2wv/L) 0[ y - 1 



x + l 



,2y-l 



(15) 



which is the Frobenius-Perron equation for the dyadic multibaker map based on quasiperiodic boundary conditions 
p4| . The respective Frobenius-Perron operator has been analyzed in detail elsewhere PJ,E5[ . Its decay rates are 



S mv {k) = In Xmu(k) 



rn In 2 + In cos k 



2-kv 



with to = 0, 1, 2, 3, v = 0, 1, 2, L — 1, and with a degeneracy of (to + 1). The eigenvectors 

(Qi D) ipmi>){x,y; k) = e Smi,(fe) ip m v(x, y;k) , 

and some root vectors have been constructed in Rcfs. [EMEa in terms of the cumulative functions 



dx' \ dy ip 0v (x' ,y'\ k) 



(16) 



(17) 



(18) 



For small enough wavenumbers k, these cumulative functions are continuous functions which are products of a mono- 
mial in x with a nondifferentiable de Rham function in y. Accordingly, the eigenvectors ipo^ are complex singular 
measures for small enough k. 

We observe that the decay rate with m — and v — vanishes quadratically as k — > in agreement with the 
expected diffusive behavior, 



In cos k 



2 



12 



(19) 



which shows that the diffusion coefficient is D = 1/2. 

The nonequilibrium steady states corresponding to a uniform density gradient g for the mean density of Eq. (^) have 
also been constructed |2(| . It has been shown that the nonequilibrium steady state of the infinite system corresponds 
to a singular invariant measure represented by a continuous but nondifferentiable cumulative function 



V 



{x,y,l) = dx' dy' f st . st .(x',y', I) = glxy + g x T(y) 



(20) 







where T(y) is the Takagi function obtained as a solution of the iteration |28[ 

i T(2y) + y , < y < \ , 



T(y) 



V 



1 , 



5<y<i- 



(21) 



The Takagi function is nondifferentiable because its formal derivative is infinite almost everywhere. It is given by a 
Lebowitz-McLennan type of formula E9| 



dT 
dy 



(22) 
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where M(y) = (2y) (modulo 1), and j(y) = ±1 if y < 1/2 or y > 1/2, respectively, is the jump -velocity. The singular 
character of the diffusive steady state turns out to be a general feature in finite-dimensional deterministic chaos of 
large spatial extension, as shown elsewhere |^2|,^6| . 

Moreover, this singular character of the steady state measure plays a fundamental role in the explanation of the 
entropy production of irreversible thermodynamics. In particular, the expected entropy production can be derived from 
the Takagi function in the case of the multibaker map |5(J . The presence of this singularity solves the famous paradox 
of the constancy of the Gibbs entropy, which can be set up when we do not recognize that the out-of-equilibrium 
invariant measure is very different from the equilibrium one on fine scales in phase space. The out-of-equilibrium 
invariant measure becomes singular if the nonequilibrium constraints are imposed at distances larger than several 
mean free paths. This is probably a paradoxical aspect of the local equilibrium hypothesis that a nonequilibrium 
system appears in local equilibrium on the largest scales of phase space although intrinsic correlations exist on finer 
scales which are due to the chaotic dynamics. The singular character of the invariant measure explains that there is an 
entropy production in large nonequilibrium systems where the chaotic dynamics removes the signature of determinism 
down to extremely fine scales in phase space. 

B. Diffusion coefficients in parameter-dependent models 

In this subsection, we describe an important consequence of deterministic chaos which shows up when a system 
parameter is varied. Under such circumstances, the diffusion coefficient of maps like a parameter-dependent multibaker 
exhibits a fractal structure by changing the parameter. In the same way as outlined in this subsection, we will show 
later that the same behavior appears in the reactive transport properties of our parameter-dependent multibaker. We 
thus emphasize the striking analogies between the diffusive and the reactive properties in regard to their parametric 
sensitivity. 

For this purpose, we summarize the main methods and results concerning the parametric variation of the diffusive 
dynamics in multibakers and in simple one-dimensional maps. To study parameter-dependent transport in multibaker 
models like Eq. (^) such one-dimensional maps are crucial, because they govern the dynamics of multibakers projected 
onto the x-axis. We therefore start with a brief review on parameter-dependent diffusion in one-dimensional maps. 
We then show how all the methods and results obtained for one-dimensional maps carry over to a two-dimensional 
parameter-dependent multibaker which we introduce and discuss to the end of this subsection. 

1. Fractal forms in a Green-Kubo formula 

Chains of one-dimensional chaotic maps are the simplest dynamical systems in which deterministic diffusion can 
be studied fl3l| . One may think of them as deterministic generalizations of simple random walks on the real line, 
where the full microscopic history of the particles is taken into account. In general, the microscopic dynamics is 
affected by changes of some control parameters. Hence, in contrast to the discussion of the previous subsection where 
no parameter has been varied, we will be interested here in the resulting parameter dependence of the diffusion 
coefficient ||[}2]-[}5) . As mentioned before, the parameters to be considered may be physically related to, for instance, 
varying the density or the shape of scatterers. An example of such systems is the chain of piecewise linear maps 
depicted in Fig. 0. This map is of the general form xt+i = M a (xt), and it is periodic by satisfying the condition 
M a (x + 1) = M a (x) + 1. The slope a serves as the control parameter and is trivially related to the Lyapunov exponent 
of the map via A = In a. The parameter-dependent diffusion coefficient can be obtained by the Green-Kubo formula 
§11111 

D^^l^aix^UMl^- 1 -^)) , (23) 

where the average (•) = f Q dx g a (x)(-) has to be taken over the invariant probability density g a (x) on the unit 
interval. j a (x) gives the integer number of boxes a particle has traversed after one iteration starting at position x and 
is thus the parameter-dependent extension of the jump-velocity introduced in Eq. ( p2[ ) above. Thus, the complete 
microscopic dynamics is divided into two parts in Eq. (p3|): the intra- cell dynamics, that is, the dynamics within a 
single box, which is represented by the invariant probability density, and the inter-cell dynamics given by the sum of 
the jump-velocities, which contains the history of the particles travelling between the single boxes of the chain. 

For computing the diffusion coefficient both parts can be treated separately. The invariant probability density is ob- 
tained by solving the Frobenius-Perron equation for the map restricted to the unit interval, M a {x) = M a (x) (mod 1), 
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Qa tt+ l(x) = I dz Q a ,t( z ) S [ x ~ M a (z)] 



(24) 



To do this, Eq. ( p4| ) can be written as a matrix equation, where the Frobenius-Perron operator has been transformed 
into a transition matrix ^,^,^2|,^6) . For maps of the type considered here, exact transition matrices can be constructed 
whenever a so-called Markov partition exists. This is the case for a dense set of parameter values a on the real line. 
The invariant probability density can then be calculated either by solving the eigenvalue problem of the transition 
matrix, which in simple cases can be performed analytically, or by solving the Frobenius-Perron equation by iterating 
the transition matrices numerically pTJ . 

In Fig. |J (a) and (b), typical invariant probability densities are plotted at two values of the slope. They are step 
functions on the unit interval, where the regions of the functions being piecewise constant correspond to the single 
cells of the respective Markov partitions. For piecewise linear maps, the invariant probability densities should always 
be step functions, although for arbitrary parameter value they may consist of infinitely many steps ]32] , |38| , |39| . The 
sum of jump- velocities, as the second ingredient of the Green-Kubo formula Eq. (p3|), gives the integer displacement 
of a particle after t iterations starting at initial position x. Since the system is chaotic, this function is highly irregular 
in x. To deal with this quantity, it is more convenient to define functions T a {x) via 



dTa 
dx 



oo 



-S S > J a [Ml{x)] 



t=0 



(25) 



which is a parametric extension of Eq. (p2|). The functions T a (x) now give the integral of the displacement of particles 
which start in a certain subinterval, and they behave much more regular in x than the sums of jumps. Employing 
T a (x) — linii^oo T a .i(x), it can be shown that these functions are obtained in terms of the recursion relation 



T a , t (x) = t a (x) + - T t-i 
a 



M a (x) 



(26) 



with t a (x) being determined by dt a /dx = j a {x) and by requiring that T a (0) = T a (l) = 0. T a (x) can be computed by 
iterating Eq. (|2^) numerically. For two special values of the slope the results have been plotted in Fig. ^ (c) and (d). 
The functions T a (x) are self-similar on the unit interval and scale with the slope a. For a = 2, Eq. ( |26|) appears as a 
special case of Eq. (pl|). Therefore, functions like T a (x) may be denoted as generalized Takagi functions. 

The numerically exact result for the parameter-dependent diffusion coefficient is shown in Fig. |5| for 2 < a < 8. 
Naively, one may have expected that D(a) increases monotonically by increasing the slope. But this is only the 
case on a sufficiently coarse grained scale, where D(a) can in fact be qualitatively matched to the results of two 
simple random walk models |33|. On a fine scale, however, D(a) shows a complicated structure with different regions 
exhibiting different kinds of self-similarity. A numerical estimation shows that the curve has a fractal dimension which 
is very close to, but greater than one. 

This highly irregular behavior of D(a) is caused by correlations of increasingly higher order in the microscopic 
dynamics of the map. For instance, in the initial region 2 < a < 3, which has been magnified in Fig. [H| the fine 
structure can be physically explained by relating local extrema on the curve to characteristics of the microscopic 
scattering process in one box as it changes with the parameter a |]J^| : if stronger backs cattering sets in by making a 
larger, the curve exhibits a local maximum, if stronger forward scattering occurrs, it goes through a local minimum. 

More generally, the fractal character of D{a) can be understood by analyzing the Green-Kubo Eq. ( p3] ) ^2 35 . Two 
basic components in the formula are responsible for the fractal character of the curve: On the one hand, the diffusion 
coefficient is given in terms of sums of jumps, which, according to Eq. (p5|), are related to fractal generalized Takagi 
functions as shown in Fig. ^ (a) and (b). This goes together with the jump velocity j a (x) having a discontinuity which 
varies with the parameter a and which reveals in a sense the fractal character of the generalized Takagi functions. On 
the other hand, a second source of irregularity are the stepwise discontinuities in the density of the invariant measure 
g a {x) as shown in Fig. |^ (c) and (d). The irregular behavior of the diffusion coefficient results from a combination 
of these effects, which are connected in the Green-Kubo formula via integrating the respective generalized Takagi 
functions over the respective invariant density. Thus, actually this behavior finds its origin in the non-robustness of 
the topology of the trajectories under parametric perturbations. 



2. A time-reversible area-preserving multibaker with fractal diffusion coefficients 

The same phenomenon of a fractal diffusion coefficient appears in a parameter-dependent generalization of the 
diffusive-reactive multibaker model introduced above |34|. This two-dimensional area- preserving map is sketched in 
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Fig. H Here, the two rectangles of the left and of the right half of the square are "sliding" along the upper and the 
lower horizontal channel of the periodically continued map governed by a parameter h, as shown in the figure. It 
should be noted that for h = 0.5 and shifting the coordinate system by Aa: = 0.5 the model reduces to the simple 
dyadic multibaker of Eq. (||). The dynamics of the probability density ft(x, y, I) of the full multibaker <&h{x, y, I) is 
determined by the Frobenius-Perron equation ft+i(x, y, I) = ft[$h (%, y, I)}, where ^^(x, y, I) is the inverse map. A 
projection of this two-dimensional Frobenius-Perron equation onto the unstable x-direction by integrating over the 
stable y-direction via g t (x,l) = J dy ft(x,y,l) |p],p3f shows that the dynamics of the probability density g t (x,l) is 
determined by the Frobenius-Perron equation of the simple one-dimensional map included in Fig. which is a kind of 
Bernoulli map shifted symmetrically by a height h. This one-dimensional map governs the dynamics of the multibaker 
map projected on the x-axis. By extending the system periodically, we recover a chain of one-dimensional maps of 
the type of the one shown in Fig. |3[ Concerning time-reversibility, we follow the definition that there must exist an 
involution G in phase space, GoG = 1, which reverses the direction of time via Go$oG = <fr 1 {I(J . For the special 
case of h taking multiples of 1/2 involutions G can be found which are related to a simple mirroring in phase space 
jyj. For general h, it can be shown that the system has strong time-reversible properties, although the existence of 
an involution G remains an open question 34,40], 



To compute the parameter-dependent diffusion coefficient of this multibaker we use that the projected dynamics is 
governed by a one-dimensional map, and thus we apply the same methods as outlined above. The result is shown in 
Fig. |^. The diffusion coefficient is again a non-trivial function of the parameter h and shares many characteristics of the 
curve presented in Fig. |^, for example, a certain random walk-like behavior on a coarse grained scale p3[ |. But it also 
exhibits some new features, especially that the diffusion coefficient is constant in intervals 0.5+m < h < 1+m, m £ No- 
This is due to the fact that the transition matrices corresponding to respective Markov partitions, and thus the 
respective symbolic dynamics of the map, do not change in this parameter interval. It is worth mentioning that in 
contrast to the specific model of Fig. || the invariant probability density of the projected one-dimensional map here is 
always uniform for all parameter values of h. Therefore, the only contributions to the fractality of D(h) come from 
the inter-cell dynamics as described by the Takagi functions T a (x). 

Along the same lines as above, we can also consider parametric variations of a bias in one-dimensional maps and 
in multibaker models In these systems, the deterministic dynamics appears in form of currents which 

are fractal functions of the bias, in certain parameter regions the mean current can run opposite to the bias, and 
the diffusion coefficient can be zero with non-zero current ]3^j3^ , ^2]| . In this regard, it is interesting to point out 
that drift currents which are irregularly fluctuating by varying the bias have also been observed numerically in other 
deterministic models H|. Moreover, we notice that certain biased maps can be related to so-called ratchets [E2U 



IV. REACTIVE DYNAMICS 



A. Reactive modes of the dyadic multibaker 

In this section, we turn to the study of the chemiodynamic or reactive modes of our simple dyadic model Eq. (|^) of 
diffusion-controlled reaction. Contrary to the total number of particles, Na + Nb, which is a constant of motion, the 
numbers of particles of each chemical species are not conserved. Accordingly, we should not expect that the reactive 
modes have a vanishing decay rate as k — > 0. This is in contrast to the diffusive modes which are related to the 
conserved total number of particles and for which the decay rate (|l9| ) vanishes at k = 0. 

Here, we consider the subspace defined by the difference between the particle concentrations in the multibaker, 

g(x,y,l) = f(x,y,l,A) - f(x,y,l,B) . (27) 

Thus, we employ the fact that the dynamics of the concentration difference can be decoupled from the mean density 
for this model, as has been mentioned before, compare to Eq. (^|). 

With quasiperiodic boundary conditions, the difference of chemical concentration evolves in time according to the 
reactive evolution operator 
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Rk= < 



g t+ i(x,y,0) = e(±-yj g t (§,2y,l) + e- ikL o(y-\) g t (^,2y - 1,L - l) 
g t +i(x,y,l) = e(\-y) <? t (f,2y,2) - d{y - §) 9t (£±I, 2y - 1, 0) , 
<7t+iOr,y,2) = ^(l-y) fft(f,2i/,3) + 5t (^,2y- 1,1) , 



(28) 



St (2±l,2y-l,L-2) 



.g t+1 (x,y,L-2) = fl(§-tf) gt(j,2y,L — l) + - |) 5t (2±i, 2y - 1, L - 3) 
j t+1 (x,y,L-l) = -e^ L e(\~y) 5t (f,2y,0) + d{y-\) 

Our goal is here to obtain the eigenvalues and eigenstates of this reactive evolution operator 

with x(^) — exp[s(A;)]. We define the cumulative functions 

G t (x,y,l) = I dx' dy' g t (x',y',l) , 
Jo Jo 



(29) 



(30) 



which obey a set of equations which can be derived from Eq. (^8|). We suppose that the leading eigenstates are 
uniform along the unstable direction x, which is justified by the fact that the hyperbolic dynamics smoothens out any 
heterogeneities along the unstable direction, 



*(x,y,l) = V{y,l) , 



(31) 



where T>(y, I) is a Schwartz distribution. We note that the further eigenstates and root states do depend on x and 
require a more detailed analysis. The cumulative functions of the leading eigenstates are thus 



G, 



cigcnstatc 



(x,y,l) = x C(y,l) with C(y,l) 



dy' V(y',l) 



(32) 



Replacing g t (x,y,l) by T>(y,l) and gt+i(x,y, I) by x£>(y, /) m Eq. fl28|) and integrating over the interval [0,y], we 
obtain the following iterative equations for the new functions C(y, I) 



C(y,0) 

C(y,l) 
C(y,2) 

C(y,L- 
C(y,L- 



& C{2yA) 



y < 1/2 , 



% [C(l, 1) + exp(-ifcL)C(2y - 1, L - 1)] , y > 1/2 , 
±C(2y,2), y<l/2, 
i [C(l,2)-C(2y-1,0)] , y>l/2, 
iC(2j/,3), 2/ < 1/2 , 
^ [C(l,3)+C(2y-l,l)] , y>l/2, 



(33) 



^C(2y,L-l), y<l/2, 
i [C(l,L-l) + C(2y-l,L-3)] 

-Sf^^O), 1/ < 1/2 , 



y > 1/2 



i [-exp(ifcL)C(l,0) + C(2y-l,L-2)] , y > 1/2 . 
The eigenvalue can be obtained by setting y = 1 in Eq. (|33|), which leads to the eigenvalue equation 

/ C(1,0) \ 
C(l,l) 
C(l,2) 
C(l,3) 

C(l,i-2) 



/ 


-2 X 


1 








■ 





exp(-ikL) \ 




-1 


-2 X 


1 





• 













1 


-2x 


1 


• 
















1 


-2X ■ 


• 






















• 1 


-2 X 


1 


V- 


exp(ifcL) 











■ 


1 


-2 X J 



(34) 
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The characteristic determinant has been calculated for several values of the distance L between the reactive sites, 



L = 


3 : 


4 X 3 


f x + cos(3fc) 


= o , 




L = 


4 : 


8x 4 


- 2 + 2cos(4fc) 


= o , 




L = 


5 : 


16 X 5 


- 4 X 3 - 3 X + 


cos(5fc) 


= o , 


L = 


6 : 


32 X 6 


- 16x 4 - 6x 2 


+ 1 + 


cos(6fc) = , 



(35) 
(36) 
(37) 
(38) 



The corresponding dispersion relations of the reactive modes are depicted in Fig. || together with those of the diffusive 
modes. Fig. || shows that the slowest decay rate which gives the reaction rate appears at k = for L odd and at 
k = ±ir/L for L even. The cumulative functions {C(y, 0}f=o °f the eigenstate corresponding to the reaction rate at 
k = are depicted in Fig. ^| for the model with L = 3 by solving Eq. ( p3|) iteratively. Near its maximum values, the 
dispersion relation behaves quadratically like 

Lodd: s (r) (£;,£) = - k(L) - D^(L)k 2 + C(fc 4 ) at fc = 0, (39) 
Leven: s {l) (k,L) = - k(L) - D {l \L) (k^Tr/lf + O Uk T ^/L)*] at k = ±ir/L. (40) 

An analytical calculation of the reaction rate h(L) and a numerical calculation of the reactive diffusion coefficient 
.E>( r ) (L) versus L reveal that 

k{L) = - In cos^ = ^ + 0(L- 4 ) , (41) 
~ i . (42) 



The reaction rate thus behaves as expected for diffusion-controlled reaction in one dimension, compare to Eq. (|1C 
These results, combined with the results for the diffusive modes, show that, on macroscopic scales, the coarse-grained 
density and the concentration difference 

p(l) = f dx f dy f(x,y,l) , (43) 
Jo Jo 

a(l) = [ dx [ dyg(x,y,l) = G(l, 1,1) , (44) 
Jo Jo 

behave like 

diffusive mode : — D , (45) 

ot oV 

reactive mode : — ~ —r — k a (L odd) . (46) 

ot ol z 

Corrections with higher-order spatial derivatives could also be taken into account in the dynamics of the reactive 
mode. For a model with L odd, this behavior corresponds to a macroscopic diffusion-reaction system with 



dp A _ D + D (r) d 2 pa D - £)« d 2 p B k 

~df ~ 2 W + 2 d! 2 2 

dp B _ D-DW d 2 PA D + D^ d 2 p B k 

~df ~ 2 W + 2 W + 2 



(Pa -Pb) , (47) 
(Pa -Pb), (48) 



where pa = p + (t/2 and pb = p — cr/2. According to these macroscopic equations, the diffusion coefficient of each 
species is Da = D b = {D + D^)/2, the cross-diffusion coefficient is Dab = Dba = (D — I)W)/2, while the reaction 
rate of Eq. (g) is given by the logarithm of the absolute value of the leading eigenvalue of the reactive evolution 
operator as 

k = | = - \ In cos ^ = - l - In \ x (k = 0)| . (49) 
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We remark that, according to the microscopic analysis, the macroscopic equations of a diffusion-reaction system do not 
necessarily follow the simple assumption often carried out that the cross-diffusion coefficients vanish, Dab = Dba = 0. 
This particular case is only recovered if D = , which is not fulfilled here. The origin of this difference holds in the 
fact that the diffusion coefficient Da associated with the state A of a particle is in general different from the diffusion 
coefficient of the particle itself which may be in two possible states A or B. In this regard, the cross-diffusion appears 
of importance in reacting systems. 

Besides, the models with L even follow more complicated diffusion-reaction equations where the reactive diffusion 
coefficient has a different status because it is associated with a nonvanishing wavenumber k = ±n/L. Nevertheless, 
the part of the diffusion-reaction process which is responsible for the reactive exponential decay is confirmed by the 
microscopic analysis. 



B. Reaction rates in the parameter-dependent multibaker 



We now discuss the parameter-dependent reactive multibaker by taking the shift parameter h into account, as it 
has already been done for the purely diffusive case (see Fig. ^|). Thus, in addition to the integer periodicity L of the 
reaction cells of the multibaker the reaction rate k will also depend on h. One may then raise the question how the 
reaction rate n(h, L) changes with respect to varying h for fixed L. Moreover, we will give some illustrative features 
of the time-dependent dynamics of the reaction process for typical h parameters. 

Analogously to the previous subsection, we start with the difference of chemical concentrations g(x, y, I) as defined 
in Eq. (|27[) . We again use the property that parallel to the x-axis the two-dimensional reactive multibaker can be 
projected onto a one-dimensional map, as has been pointed out before (see Fig. ^|). The time evolution of the projected 
reactive part (,t(x,l) = J dy g(x, y, I) of the multibaker is then determined by the reactive evolution equation of a 
respective one-dimensional map, 

C t (x,o = i? (1) (/*,£R t -i(z,0 ■ (50) 

Here, R^(h,L) represents the one-dimensional reactive evolution operator, and Q(x,l) = Qt(x,l,A) — g t (x,l,B) is 
the difference between A and -B-particle densities in the corresponding one-dimensional map. As has been done for 
the purely diffusive case, we again write this equation as a matrix equation, where instead of R^'Qi, L) a topological 
transition matrix T(h, L) acts onto a particle density vector t. The matrix T(h, L) is structured such that in case of 
reactive scattering centers the elements in the corresponding columns of the matrix have a negative sign, and thus a 
particle changes color by leaving a reaction cell. Otherwise, the matrix is the same as discussed for the diffusive case. 

We first discuss some details of the time evolution of the reactive modes. By integrating over C*(x, I) or its respective 
vector representation we obtain the difference between the total number of A and _B-particles at discrete time t which 



is £ t = Yli I dx ( t (x,l). From the corresponding phenomenological time-continuous reaction equation Eq. (46) one 
would expect that for the reactive multibaker £ t decays exponentially after a suitable coarse graining according to 
£t = £o ex p(~ Ht). If this is the case, we can define the reaction rate of the reactive multibaker in analogy to the 
phenomenological equation as K(h,L) = k/2. 



To compute n(h,L) according to this definition, we solve the matrix formulation of Eq. (50) by iterating the 
transition matrices T(h,L) numerically. As an initial particle density we choose C^o(x,l) to be uniform in one reactive 
cell of the multibaker, which corresponds to having only A particles in this cell with the number of B particles being 
locally zero, and we make the chain long enough such that the evolving density is not affected by boundary conditions. 
Fig. (a) and (b) give two typical examples of (t(x, I) for certain parameter values of h after t — 40 iterations. They 
show how the "perturbation" (,q(x, I), which is a local initial deviation from the equilibrium state Ct( x i = (t — ■> oo 
spreads out along the rr-axis by exhibiting a rather complex fine structure with oscillations around zero. Fig. R 
(c) and (d) contain half-logarithmic plots of |£ t | with respect to the discrete time t. These plots reflect a different 
dynamical behavior of £ f for different magnitudes of the reaction rate. For n(h, L) close to zero, see the upper two 



curves in Fig. 10 (c), £t decays apparently non-exponentially for small times t. Only for larger times it eventually 
reaches exponential decay. Thus, the system shows that it is close to states of the h parameter where it is non-reacting. 
For the lowest curve in Fig. |To| (c), which corresponds to an intermediate reaction rate, £t provides initially strong 
periodic fluctuations. They are partly due to the complex deterministic dynamics of the reactive baker in one cell of 
the chain, as has already been observed and explained for a one-dimensional purely diffusive case Apart from 
such strong periodic oscillations on a fine scale, in Fig. ^ (d) £t exhibits an interesting crossover between a fast decay 
for smaller times and a slower decay for larger times, where again it approaches exponential behavior. This may 
reflect the fact that for larger reaction rates n(h, L) ^> the reaction is controlled by the diffusive dynamics. These 
features of £t should be compared to the characteristics of the respective probability densities Ct(%> in the figure. 
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To obtain quantitative values for the reaction rate, Eq. ( |50| ) has been analyzed by solving the eigenvalue problem 
of the corresponding transition matrix T(/i, L), analogously to what has been done in the previous subsection for the 
dyadic reactive multibaker. In general, the spectra of T(h,L) are extremely complicated [Q. However, as has been 
argued in the previous subsection for special cases, and supported by our observation of long-time exponential decay 
of £t for the general case, we expect that, in the limit of infinite time, the reaction in the multibaker will always 
be governed by the slowest eigenmodes and their respective eigenvalues. This motivates to define the parameter- 
dependent reaction rate kQi, L) via the maximum of the absolute value of the eigenvalues of T(/i, L), 

K(h,L) = ~]n\xmUh,L)\ , (51) 



analogously to Eq. (49). Numerically, we find that for large regions of the h parameter a certain fundamental domain 
Lf of the multibaker is sufficient to obtain the correct leading eigenvalue Xmax(h, L). This domain must always 
include multiples of two reactive centers, and its length is defined by the number Lp = 2L Int(l + h) of cells of 
the multibaker. In these regions, solutions for the eigenvalue problem of T(h,L) defined on the domain Lp lead 
to a maximum eigenvalue Xmax{h,L) as obtained by solving the corresponding eigenvalue problem for longer and 
longer chain lengths mL — > oo , m 6 N. However, especially for small L and large kQi, L) this fundamental domain 
only provides an approximation to the exact results which are then obtained by making the chain length mL large 
enough such that the error in Xma,x(h, L) with respect to mL is sufficiently small. Fig. O shows some typical largest 
eigenmodes tj){x,l) on the fundamental domain in cases where it gives the correct corresponding largest eigenvalue 
K,(h,L). For large reaction rates, to a certain respect the largest eigenmodes behave like sine functions, see Fig. |ll| 
(a), whereas for smaller reaction rates, the largest eigenmodes approach two-periodic step- like functions as shown in 
Fig. [ll] (b). Fig. [ll] (c) depicts the largest eigenmodes for a parameter value of h which is just at the borderline of a 
non-reacting h region, but where the system is nevertheless already highly reactive. Here, the eigenmodes appear to 
be especially complicated p8[ . 

In Fig. [jj] (a) the reaction rate n(h, L) as defined via Eq. ( pl| ) has been computed for a series of reaction center 
periodicities L. For h — 0.25 there is no reaction rate in the system. In this case, the iteration method confirms 
that the difference in the number of particles £t oscillates periodically around zero instead of decaying exponentially. 
An analysis of the eigenvalue spectra of the corresponding transition matrices reveals that at this h parameter the 
respective reactive multibakers are not ergodic p£|. For all other h of the figure the reaction rate is well-defined and 
shows a complicated structure. By increasing L the reaction rate kQi, L) decreases almost everywhere, as one can 
expect intuitively, except in certain small parameter regions of h. Fig. |l2| (b) gives the full result for n{h, L) at L = 2. 
This structure repeats itself with a periodicity of 2m < h < 2 + 2m , m £ No- In certain intervals of h the four peaks 
depicted in the figure are very similar, or even identical, however, there does not appear to be a simple scaling law by 
which the full peaks can be mapped onto each other. The plateau regions with zero reaction rate correspond to the 
respective regions observed in Fig. ^ for the diffusion coefficient of the system. They share the same characteristics 
as discussed above for the singular case of h = 0.25, except that at h = 1 the system is ergodic, but not mixing. 
Topologically, these regions are of the same origin as explained for the diffusive case. In Fig. (c) the reaction rate 
has been computed for L — 3. In contrast to the two-periodic case the change from a non-reactive region to a reactive 
region occurs for L = 3 apparently discontinuously in the reaction rate by varying h around 1. This corresponds to 
the system at h = 1 and L — 3 being mixing, whereas for the same h and L = 2 it was ergodic, but not mixing. 
We note that for L = 4 there are even two of such discontinuous transitions. The detailed irregular structure of 
the curves, as well as the phase-transition like behavior shown in Fig. [l^ (c), can be understood in more detail by 
analyzing the eigenvalue spectra of the reactive evolution operator and how they change under parameter variation, 
as will be discussed elsewhere Q . 

Apart from varying h, other parameter dependencies can be studied in this reactive multibaker as well. For example, 
the distance between the single reaction centers could be changed by allowing L to be continuous, the positions of 
the reaction centers could be shifted by keeping L fixed, and the size of the reaction centers could be increased or 
decreased. In all these cases we expect non-trivial parameter dependences to be typical which are similar to the one 
depicted in Fig. O @. 



V. DISCUSSION AND CONCLUSIONS 



In this paper, we have analyzed simple deterministic models of diffusion-controlled reaction. The models fulfill the 
chaotic hypothesis mentioned in the Introduction, which allows us a much sharper analysis of the phenomcnological 
foundation of diffusion-reaction processes than with the old stochastic assumption. In this regard, for the simplest 
model we have been able to derive the exact dispersion relations not only of the diffusive modes but also of the reactive 
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modes. The reactive modes indeed decay exponentially, as it should be for nonconserved quantities. The reaction rate 
behaves like k ~ D/L 2 , as expected for diffusion-controlled reactions in one dimension. The reaction rate introduces 
a new time scale T react ~ k^ 1 ~ L 2 /D, with respect to purely diffusive systems. This reactive time depends on 
the concentration ~ l/L of catalysts and, therefore, it takes a fixed value in a given diffusion-reaction system. This 
time scale of reaction is intermediate between the short time scale of the Lyapunov instability T c h a0 s ~ A^ x and 
the long time scale of diffusion T^g ~ C 2 /D. An important difference between diffusion and reaction is that the 
evolution operator is positive in the diffusive subspace and is thus of Frobenius-Perron type, although the evolution 
operator (|28|) has both signs in the reactive subspace. By analyzing this deterministic evolution operator, we have 
found that the spatial dynamics of the reactive modes appears significantly different from the standard assumption 
of macroscopic diffusion-reaction models due to the particular importance of cross diffusion. On the other hand, the 
eigenstates associated with the reactive and diffusive modes are expressed as singular Schwartz distributions, also in 
contrast with the phenomenological models which have always suggested that the eigenstates are regular functions. 
In a sense, we may say that this singular character of the exact eigenmodes renders their relaxation compatible with 
the deterministic dynamics of the particles, in full respect to the mechanical Liouville theorem of volume preservation. 
This result previously observed for diffusion is here shown to hold also for reaction. 

Moreover, the macroscopic transport coefficients such as the diffusion coefficient and the reaction rate both turn 
out to exhibit a highly irregular behavior as a function of a control parameter of the system. In this way, we 
have pursued the early work by Elskens and Kapral who studied the isomerization rate for the simplest values of 
their model parameter |f[gj| . The irregular behavior has its origin in the topological instability of the trajectories 
in phase space and also in the singular character as represented by Takagi functions. It appears as a fundamental 
result because the topological instability of phase-space dynamics is a common feature to many systems including the 
Lorentz gas, the hard-sphere gas, and nonhyperbolic systems which are all non-robust under parametric perturbations. 
However, we may expect that the transport coefficients would have an irregular behavior only if the systems fulfill 
certain additional criteria, as being spatially periodic, being sufficiently low-dimensional, and being such that particle- 
particle interactions are not of main importance. Physical systems of this kind could - to a certain extent - already 
be realized experimentally in form of so-called antidot lattices . On the other hand, the detailed fractal character 
of parameter-dependent transport coefficients may disappear by increasing the number of degrees of freedom, or by 
including strong stochastic perturbations . We have moreover observed that the reaction rate is not only a highly 
irregular function of the parameter but also presents discontinuities which are reminiscent of nonequilibrium phase 
transitions. The parametric sensitivity seems thus enhanced at the level of the reactive properties as compared with 
the diffusive ones. 
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FIG. 1. Examples of reactive Lorentz-type models: (a) on a regular lattice; (b) on a random lattice. 
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FIG. 2. Geometry of a reactive multibaker or bakery map. 
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FIG. 4. (a), (b) Invariant probability densities g a {x) on the unit interval for the map of Fig. |§] modulo 1. The slope is 
a ~ 2.5004 for (a) and a ~ 3.49997 for (b). (c), (d) Generalized Takagi functions T a (x) for the same map at a = 3 in (c) and 
at a = 4 in (d). 
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FIG. 5. Parameter-dependent diffusion coefficient D{a) for the map of Fig. ^| and blowup of the initial region. The main 
jraph consists of 7,908 single data points, the magnification of 979. In both cases errorbars are too small to be visible. 




FIG. 6. Dynamics of one cell of an area-preserving time-reversible multibaker with a non-trivial parameter-dependence h. 
Projection of the dynamics onto the horizontal axis reduces the system to the symmetric one-dimensional piecewise linear map 
shown in the figure to the left which for h = is the Bernoulli shift. 
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FIG. 7. Parameter-dependent diffusion coefficient D(h) for the multibaker of Fig. ^ and blowup of the initial region. The 
main graph consists of 638 data points, the magnification of 514. In both cases the single points have been connected with 
lines, errorbars are too small to be visible. 




FIG. 8. Dispersion relations of the diffusive modes (solid lines) and of the reactive modes (dotted lines) for the dyadic 
reactive multibaker with (a) L = 3; (b) L = 4; (c) L — 7; (d) L — 8. (note: poor quality output; the originals can be obtained 
upon request) 
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FIG. 10. (a), (b) Densities corresponding to the difference in the number of A and B-particles per total number of particles 
after t — 40 iterations of the map. In (a) (L = 2, h ~ 0.2429) the reaction rate is close to zero, whereas in (b) (L — 3, h ~ 0.1496) 
it is locally maximal in h. (c), (d) Half- logarithmic plots of the total difference £t in the number of A and B-particles of the 
system as it varies in time t. In (c) the corresponding reaction rates n(h, L) are close to zero for the upper two curves (k < 0.006), 
for the lower curve the reaction rate has an intermediate value (k ~ 0.02), whereas in (d) it is locally maximal in h (k ~ 0.05). 
The parameters for the upper curve in (c) correspond to (a), the curve in the middle is at L = 3, h ~ 0.247, and the lowest one 
is at L — 3, h ~ 0.4472. Case (d) corresponds to (b). 
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FIG. 11. Examples of largest eigenmodes t/j(x,l) for the reactive multibaker corresponding to the largest eigenvalue 
Xmax(h,L) parallel to the x axis in the fundamental domain L F as described in the text. For (a) (L = 2, h ~ 0.1496) 
the reaction rate n(h,L) is very large (k ~ 0.12), for (b) (L = 2, h ~ 0.4947) it is very small (k ~ 0.003). In both cases, 
there exist only two real largest eigenmodes where the second ones are shifted by a phase. In (c) (L — 6, h = 1) both largest 
eigenmodes have been plotted (thick line for the one and thin line for the other, respectively). 
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